function yhat = SSVM_LossAumentedInference(w, x, y)
  %fprintf('constraintCB x=[%d %d],y=%d\n',x,y);
% slack resaling: argmax_y delta(yi, y) (1 + <psi(x,y), w> - <psi(x,yi), w>)
% margin rescaling: argmax_y delta(yi, y) + <psi(x,y), w>
  if dot(y*x, w) > 1, yhat = y ; else yhat = -y ; end
  
end

